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1. INTRODUCTION 


The purpose of this paper is to explore the use of lunar surface 
photography in order to achieve the photogrammetric transfer of available 
selenographic coordinates from future lunar landing sites to neighboring, 
photoidentifiable features. It can be implied from the procedures developed 
that overhead photography, were it available, could be utilized and would 
provide a material strengthening of the total solution. By the methodic 
selection of features and confirmation that they can in reality be identified 
from orbital photography, a modest selenodetic control system can be 
expanded into a net that could ultimately control all future, manned or 


unmanned, orbital photographic missions, 





2. ABSTRACT 


For centuries man has scrutinized the moon in one manner or 
another and postulated theories concerning its size, shape, origin, and 
other general characteristics. With the passage of time and the improve- 
ment of equipment and observation techniques the desire for more explicit 
information concerning earth's nearest celestial neighbor has become acute. 
In fact, as the moment approached when man would actually set foot on the 
lunar surface, the need for such information became vital. The following 
historical review briefly outlines man's effort to improve his knowledge 
in one of the pertinent regions of selenodesy — selenodetic control. 

The remainder of this paper explores a method of improving the 
existing selenodetic control by employing available lunar surface photography 
supplemented by that obtained from lunar orbit. Following the results of 
this experiment an ideal model is submitted. The unknowns associated with 
this model are perturbed within realistic limits by a random number gener- 
ation program. This provides a theoretical indication of the accuracy that 
could be anticipated assuming there is reasonable adherence to the suggested 
procedures. 

Finally, conclusions are drawn and reasonable recommendations are 
offered to improve selenodetic control by the photogrammetric transfer of 
known or assumed,local or astronomic coordinates of a lunar landing site to 


neighboring features that may be photoidentified from orbital photographs. 





3. HISTORICAL REVIEW 


For decades in the past to the present day the task of surveying the 
moon has engaged the efforts of many astronomers. In early 1959 the 
launching of LUNIK I by Russia, and subsequently, POINEER IV by the 
United States: '"... opened the first modern, post telescopic phase of lunar 
exploration. ''or, at least, introduced a tantalizing new dimension [30]. 

During the some seventy years prior to the launching of the first 
lunar space probes the establishment of selenodetic control was founded on 
direct astronomic angular observations and indirect angular observations 
through astronomic photography. Essentially, it was based on heliometric 
observations which consist of measurements of position angles and angular 
distances between a reference point on the lunar surface (Mosting A is the 
fundamental point) and the lunar limb. Observations at mean libration 
permit a best-fit circle of the lunar disc to be established. The center of 
this circle is defined as the projection of the origin of the coordinate 
system (the dynamic or mass center of the moon) upon the lunar surface 
and its radius to be the mean radius of the moon. Thus, the center of 
figure is equated to the center of mass and in the adjustment of the helio- 
metric observations this injects the so-called center of figure bias. The 
adjustment provides corrected values of physical lunar libration parameters 
and the coordinates of the reference point as well as the mean radius of the 
moon [26]. 

The heliometer was first developed by Bouger in 1748 and later 
modified by Dollond. It consists of a refractor telescope with two semi-lenses 
which may form a single, superimposed image of two object points at the 
principle focus. The angular distance between the two object points formed 
on the focus is equal to the-distance between the centers of the semi-lenses 
when one slides parallel to a line of section upon the other [24 |. It was used 


to measure the diameter of the moon at the end of the 18th century byLalande 
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and by Bessel in 1839 to investigate lunar physical librations. It was 
Bessel that developed the procedures for measurement that remain basically 
intact today. 

Heliometric observations are limited by the resolving power im- 
posed by their relatively small aperatures (4-7 inches). The Rayleigh 


criterion: 


i- A 
G6=1.22 D 


§@ = minimum angle resolved in minutes 
dX = wavelength of light 


D = diameter of objective lens 


theoretically indicates that a six inch aperature provides a minimum resolu- 
tion of 0.75 arc seconds or well over a kilometer on the moon's surface. A 
further limitation is based on atmospheric refractivity [20 }. 

Selenodetic control systems derived from earth-based lunar photo- 
graphy generally rely heavily on heliometric observations. The reduction 
of these observations provide the libration parameters (f, I) and the coor- 
dinates of fundamental points. These reference points provide the orientation 
and scale of the photographs from which the plate constants are determined. 
| A German astronomer, Franz, established the original eight funda- 
mental points in the early 1900's. Through the use of five plates from Lick 
Observatory he expanded these to a system of 150 points. By 1958, an 
Austrian astronomer, Schrutka-Rechtenstamm, published a revision of the 
moon libration theory and a recomputation of Franz's 150 points. This 
system is considered the best available and has served as the basis for later, 
more densified systems [26], Yet the S-R system and others comparable to it 
reflect the inaccuracies inherent in the original heliometric observations as 
well as the additional inaccuracies associated with the earth-based photographic 
process. 


Two American government agencies have undertaken densification of 





lunar control. The Army Map Service (now, Army Topographic Command) 
published AMS-64 consisting of 256 points. This agency utilized the funda- 
mental points from the IAU Cataloque of Blagg and Muller and plates from 
the Lick Observatory [8]. In 1966, AMS published the GROUP NASA system 
of 484 points utilizing control points determined by Saunder, Franz, and 
Konig [18]. The Aeronautical Chart and Information Center of the U.S. Air 
Force published another independent system of 196 points in 1965. ACIC 
selected Control from the S-R system and plates from the Pic du Midi 
Observatory in France andthe U.S, Navy Astrometric Reflector in Arizona 
[23]. There were large differences betwcen the systems of the two agencies in 
planimetry (several kilometers) and height. This was emphasized during the 
RANGER probes to the moon when elevation differences of approximately 2.5 
kilometers between the AMS/ACIC systems and the trajectory computations 
were noted. Nevertheless, the systems were combined to form the Selenodetic 
Control System, DOD-66, of 734 points [26][19]. 

Two modern photographic methods are independent of control establish- 
ed through heliometric observations and appear to be rather promising. The 
Lunar and Planetary Laboratory at the University of Arizona employs a pro- 
cedure using star trailed photography that was designed by Arthur [26]. Per- 
haps more significant is a procedure contributed by Kopal of the University of 
Manchester. Moutsoulas describes it as photographing a stellar field that 
is at the same declination and hour angle that the moon will attain at a later 
time. When the moon reaches the proper position, the plate is reexposed. 
Providing no excessive temperature changes take place during the period the 
telescope is stationary, the star field provides the plate orientation and scale; 
and the constants can be used for reduction of points on the lunar surface [24 ]}. 
Kopal states that the achieved accuracy is sufficient to determine the physical 
librations of the moon [22]. 

Extensive, extraterrestrial photography was inaugurated with the 


launching of the Lunar Orbiter Satellites during the period August 1966 and 





August 1967. The mission of the first three Orbiters was primarily 
designed for the selection of primary and secondary landing sites for sub- 
sequent Apollo missions. Orbiter IV and V were tasked to perform a broad, 
systematic survey of scientifically interesting features on the lunar surface. 

All Orbiter photographic subsystems contained a medium resolution 
lens (focal length 80 mm) and a high resolution lens (focal length 610 mm). 
Neither was of photogrammetric quality. Calibration of the system, in 
general, included determination of the calibrated focal length, radial and 
tangential distortion, the principle point of autocollimation and the camera 
format reference system with respect to sawtooth fiducials and a preexposed 
reseau system on the film (Lunar Orbiter I lacked these reseau marks). Ad- 
ditional calibration was required to establish the effect of an image motion 
compensation system. 

In operation, the film would be clamped to the platen, and the platen 
would move in proportion to ground speed while the shutter was open. The 
film was then processed by a BIMAT system which developed, fixed, and 
dried it. The negative was then scanned by a line scan tube in small increments 
(2.67 mm). This signal was electronically processed for transmission to 
earth via the spacecrafts' telemetry subsystem as a composite video signal. 
The ground reconstruction electronics system received the video signal and 
fed it to a kinescope tube from which it was copied on 35 mm television record- 
ing film. A reassembly printer utilized this record to orient and project the 
framelets on aerographic duplicating film to produce the finished product. 

The photography collected from this series eliminated several signif- 
icant limitations attached to earth based photography; namely, the distortions 
associated with atmospheric refractivity and insufficient scale for effective 
resolution. Further, it provided a greatly improved geometry. However, 
other disadvantages inherent in the total system design requirements introduced 
distortions into the photography and uncertainties into the reduction procedures. 
Broadly, the distortions were associated with on board photographic processing, 
space transmission of the video signal, and ultimate reconstruction of the photo. 
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Reduction uncertainties included the film distortion, but additionally, was 
largely dependent upon photo support data which defined spacecraft location 
and attitude at time of exposure. These were functions of the orbit determin- 
ation program with its associated uncertainties. 

Nevertheless, despite the fundamental inaccuracies, ACIC evaluated 
the feasibility of establishing a lunar geodetic system from Lunar Orbiter 
photography and arrived at positive conclusions[3]. One result was, ae 
Positional Reference System of Lunar Features Determined From Lunar 
Orbiter Photography. Although the original feasibility study encompassed 
only the Lunar Orbiter IV Mission with its polar orbit and extensive coverage, 
it was found that the medium resolution photography was of particularly poor 
quality in detail except near the terminator. The remainder was either 
highly over or under exposed. All photography possessed significant errors 
in timing, exposure orientation, and spacecraft positioning/3]. As a result, 
photography from all Lunar Orbiter missions was utilized in order to achieve 
the desired coverage. However, Lunar Orbiter I photos which lacked a pre- 
exposed reseau grid onthe film were employed only when necessary to fill in 
specified areas. The method used, broadly, for this control system is best 


described by the author: 


"The method consists of computing perspective projections [23] 
based upon the orbital data for a series of photographs that 
are linked together by common coverage. Starting on the 
nearside | of the moon], the projections were positioned 
to agree with the coordinates of features determined 
from telescope photography. |The ACIC net of 196 control 
points,{23]], The link was continued around the moon by 
extending the coordinates of common features from one 
photograph to the next. A meridional arc and an equatorial 
arc were completed and joined in the vicinity of the equator 


and the 180th meridian [27]. 
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This net produced (considering the extent of the net and lack of farside 
control) reasonable estimated accuracies of 1-5. 5-10, and 10-15 kilo- 
meters, depending on the particular area cited [27]. This was achieved 
despite the facts that control was provided only on the nearside in a coor- 
dinate system based on center of figure, and the photography was of variable 
quality with all the errors associated with its on board processing and elec- 
tronic transmission. Further, the exterior elements of camera orientation 
were determined from spacecraft telemetry with the associated orbit deter- 
mination uncertainties and a coordinate system originating at the center of 
mass. 

A current control net in the process of being established by the 
Mapping Sciences Branch of the Manned Spacecraft Center, Houston, Texas, 
is in the imminent stages of completion. This net is based on medium and 
high resolution photography acquired solely by Lunar Orbiter IV. It covers 
a rather extensive area between +20° latitude and 60° west longitude to 45° 
east longitude with the greatest concentration of control in the Apollo land- 
ing zone of +5~ latitude of the same area. Although control points from 
DOD-66 and the ACIC/AMS nets are input to the computational program, they 
are generally not used in the adjustment. They are merely compared to the 
control established by Lunar Orbiter IV and the root mean square differences 
are output in the statistical summary. Preliminary results have shown a 
bias between the two systems of approximately two kilometers, but the final 
results have yet to be published. 

All of these control systems are steps toward the fulfillment of the 
essential requirements for the development of geodetic and cartographic 
knowledge of the moon as outlined by the Falmouth conference of scientists, 
convened by the National Aeronautics and Space Administration at Falmouth, 


Massachusetts in July of 1965 [12]. Among these requirements are: 


Establish a selenodetic coordinate system... related 


to the right ascension/declination system. 





Derive a reference figure with respect to a point 
which is representative of the moon's center of 


mass. 


Establish a three-dimensional geodetic control 
system...in terms of latitude, longitude, and 


height above the chosen reference figure. 


These requisites are not only essential to the expansionof geodetic and 
cartographic knowledge of the moon, but become fundamental, base know- 
ledge for the exercise of other disciplines [12]. Photogrammetry has 
demonstrated uniquely that it provides the necessary capability to efficiently 
gather the necessary data and to process it into useful and meaningful infor- 
mation 12]. 

The following photogrammetric procedure is submitted as a modest 
contribution to the ever expanding numbers of methods designed to increase 


man's knowledge of the lunar body. 





4, EXPERIMENTATION 


4.1 Real Data Experiment 


4.1.1 Preliminary 


The purpose of this demonstration is to describe in detail the 
procedure utilized to transfer local or selenographic coordinates from an 
assumed or known location to surrounding lunar features that are identi- 
able in orbital photographs. It must be realised, however, that no lunar 
surface photography has been accomplished with this purpose in mind. As 
a result several basic assumptions are employed and various procedures 
inaugurated that would normally be unnecessary were such a mission 
assigned to personnel of the APOLLO series or follow-on series which would 


reach the lunar surface. 


4.1.2 Materials 
The following materials, equipment, and systems were used: 
A. APOLLO 12 Lunar Surface Photographs; AS12-48-7090, 7091, 
7092; Magazine X; Exposed by a 70mm Hasselblad camera 
with focal plane reseau grid. (Nominal focal length, 60mm) 
B. A.M.S., Lunar Map, Surveyor III Site; Scale; 1:2000 
(lst ed., Jan 1968) 
C. Mann Precision Comparator, Type 735 with Mann Data Logger 
D. IBM 360/75 Computer System (OSU installation) 
The photographs identified in A. above were the result of an extensive 
search of all surface photography obtained furing the surface operations of Apollo 


Missions XI and XII. They were selected with the following criteria in mind: 


A. Stereoscopic coverage 

B. Maximum base between photographs 

Cc. Simultaneous, photographic coverage of the LM, Surveyor IT; 
and other points on the lunar surface that could be identified 


from orbital photography 
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D. Exposed with a calibrated camera equipped with a focal plane, 

reseau grid. | 

These three photos fulfilled these requirements adequately with an 
average base estimated to be twenty meters; the LM and Surveyor III were 
imaged on each photo; three relatively well defined lunar features were 
imaged; and a post flight calibration was conducted on the two cameras em- 
ployed. Each camera was equipped with reseau grid at the focal plane. 
Unfortunately, it has not been ascertained which camera exposed these par- 
ticular plates [4]. However, their calibrated focal lengths of 61.547mm 
(#1016) and 61. 636mm (#1002) determined at a 22.5m focus with black and 
white film (KODAK S0267) were quite similar [5][6]. Neither camera had a 
lens distortion pattern that would require consideration except for the most 
rigorous photogrammetric procedures [5 ][6]. 

For the purpose of this demonstration the average focal length was 
used in calculations. This constituted the introduction of approximately 
+0.07% error in the focal length and a proportional amount in the computations 
associated with it. This was considered insignificant for the purpose of the 
real data experimentation. Further the reseau grid was assumed to be at 
exactly spaced internals of 10mm, (4), and radial and tangential distortions 


were neglected [17 ][5 [6]. 


4.1.3 Procedure 

Broad exposure to the many hundreds of photographs taken during 
APOLLO XII surface operations permitted the viewer to acquire a semblance of 
orientation in regard to several features on the lunar surface. This was not 
facilitated by any documentation concerning time, direction of exposure, orien- 
tation of the camera or any other details except in the most general sense. 
Nevertheless, this orientation permitted the selection of three photographs with 
the LM, the Surveyor III and three other photoidentifiable features which could 
be located on the lunar maps. Further, it was confirmed that these features 
could be seen on available orbital photography. Specifically, this was photo- 
graphy from Lunar Orbiters I and III. APOLLO XII orbital photography which 
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covercd Surveyor III Site was taken at a height of approximately 60 nautical 
miles using a lens of 80mm focal length. The comcomittant photo scale 
was nearly 1:1,400,000. This was entirely inadequate for seniaee feature 
identification within the limitations of surface acquired photography. 

The lunar maps of Site III were employed to establish the coor- 
dinates of the five points to be used. The LM was plotted on Lunar Map, 
Surveyor III Site (Scale; 1:2000) from coordinates established on Lunar 
Surface Exploration Map, LSE 7-6, Scale 1:5000, prepared by the U.S. 
Army Topographic Command, 1 November 1969. With the top, center of the 
LM arbitrarily defining the origin of a local cartesian coordinate system its 
azimuth from Surveyor III was measured on map B as 301° 30° 00/0 and 
fixed to establish orientation. Additionally, the distance between the LM and 
Surveyor was measured and fixed at 202.00 meters to establish scale. The 
local coordinates of the three other points were obtained relative to the LM. 
The heights were determined relative to the top center of the LM by inter- 
polating between the five meter supplementary contour intervals provided on 
the map. The initial locations of all points are summarized as follows (See 


Figures I, IA, and II): 


SE LENOGRAPHIC COORDS, LOCAL CARTESIAN COORDS, 


POINT LATITUDE LONGITUDE Ok ne LZ 
1 (LM) 3-11-51.6S 23-23-14.60W 1000.00 1000. 00 100. 00 
2 (SURVEYOR) 3-12-04.0S 23-22-53.6W 1172.23 894, 46 87.49 
3 (MOUND) 3-11-46.15 23-23-20,.3 W 948, 00 1045, 00 93.96 
4 (LONE ROCK) 3-11-52.9S 23-22-58.8W 1129, 23 988, 46 93.96 
5 (CRATER RK) 3-11-53.5 S 23-22-55.7W 1156, 23 982,46 91, 96 


The location of camera exposure stations provided a more difficult 
problem since there was no documentation in their regard. Therefore, esti- 
mated positions had to be determined from the photographs themselves. This 
was accomplished graphically be constructing a template based on the camera 
field of view. With a nominal focal length of 60 millimeters and usable camera 


format of 52 by 52 millimeters the angular field of view was computed to be 
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approximately 46°. There was an angular field of 9°2 between adjacent 

reseau crosses. The template was overlayed on the lunar map and adjusted 
until identifiable lunar features were in their proper angular relationship. 

When the optimum fitting of the template was achieved, the vertex defined 
approximations of the exposure station in planimetry (X,, Y,) and the central 
axis of the template defined the direction of the camera optical axis. This 
provided an estimate for the phi (©) rotation. Exposure station height (Z,) 

was again interpolated from contour intervals modified by an added 1.37 meters 
based on the assumption that the astronaut accomplished the photography stand- 
ing with the camera at mid-chest level. Estimates of the omega (w) and 

kappa (x) rotations were determined from the apparent depression angle of 

the center cross reseau and the comparison of a line of horizontal reseau marks 
with the apparent lunar horizon, respectively. A summary of the locations of 
the exposure stations and camera orientation estimates are (See Figures I, IA, 


and Ii): 


STATION SELENOGRAPHIC COORDS. LOCAL CARTESIAN COORDS. 
(PHOTO #) LATITUDE LONGITUDE Xo Yo ie 
1 3°12/1143 § 23°22'52”0 W 1186.23m 832.46m 94,09m 
7” 3.12 09.0S 23 22 49.6W 1206.23 852.96 94,24 
a 3.12 06.7S 23 22 48.6W 1214.73 871.46 95.34 
(7092) 
ORIENT ATION (DEGREES/RADIANS) 
x 2 G2 
1 3.50 /0.06109* 20.0 / 0.34907 80.0 / 1.39626* 
2 3.50 / 0.06109 42.0 / 0.73304 80.0 / 1.39626 
3 3.50 / 0.06109 60.0 / 1, 04720 80.0 / 1.39626 


* A selected average for the three photographs was imployed for the x and 
w rotations. | 
- It became apparent during the template fitting procedure that there 
existed a definite possibility of a significant discrepancy between the location 
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Figure I. (Photo AS12-48-7092) 


16 





an an ae 





plotted for the LM and the position indicated by its angular relationship with 
other features. It appeared that the actual position of the LM should be some 
25 meters to the NE of its current position. However, since no better 
information on its selenographic coordinates was available. It was consid- 
ered to be fixed with the qualification that this discrepancy would be investi- 
gated by varying the application of constraints on it and other points during 
the adjustment. 

The original intention was to measure photo coordinates on the 
Zeiss, Precision Stereocomparator, PSK, with ancillary IBM 026 card punch 
to facilitate use of the computer program COMCORDCON, This program 
converts comparator coordinates to photo coordinates by an affine transfor- 
mation, simultaneously correcting for lens distortion and film shrinkage 
(See Appendix I). Because of the malfunction of this equipment the Mann 
Precision Comparator was utilized. Unfortunately, to simplify the obser- 
vation procedure, each plate was rotated approximately 30° to prevent 
alignment of the measuring cross with the photographic reseat crosses. This 
prohibited COMCORDCON from properly identifying the four reseau marks 
associated with each point measured and correlating them to the reseau, photo- 
coordinate system. A simple, two-point transformation routine was employed 
to rotate the comparator coordinate system near enough to the reseau photo- 
coordinate system to make the data compatible to COMCORDCON. The out- 
put from COMCORDCON was then ready for input to the BLOCK TRIANGULA- 
TION computer program (See Appendix II). 

The following mean standard errors were estimated for conversion 
to the variance-covariance matrices for subsequent use in the BLOCK 


TRIANGULATION program for weighting: 


Photo coordinates; 6, =6,= 0.01 mm 
Exterior orientation; O,, = oe = 6,, = 20-0m 


eae Ona oe a meacde(tl0-) 


aA 


Oy = 0.08727 rad (5°) 


aA 


Survey coordinates; Oo; = a= 6G, = “to 0.01 m (various) 
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4,1.4 Real Data Results | 

In addition to the variance-covariance matrices postulated from the 
standard errors of the previous section, constraints on the survey coordinates 
of Points 1 and 2 and the elevation coordinate of Point 3 were imposed assum- 
ing a standard error of 0.01 meter. The results of this first adjustment 
were exceedingly poor. Subsequent adjustments consisted of input imposing 
constraints on combinations of Points 1 and 2 and variable constraints and 
relaxations on Points 3,4, and 5. These triangulations either provided only 
slightly improved results or the adjustment failed to converge at all. 

Two tendencies were manifest, particularly. Point 1, the LM, 
continually drifted to the lunar northeast or east, and there was a constant 
warping of the model most evident in the residuals on surface point elevations, 
the x rotation which was constrained to 5°, and in w which was constrained to 
10°. When the constraints on Point 1 were relaxed, the LM freely moved 
approximately 48 meters almost due east of its initially plotted position. 

The warping appeared to subside to some extent, but further variations of the 
weight matrices were required to reduce the residuals on survey elevations 
and the rotations associated with the elements of exterior orientation to any 
degree of realism or consistency. 

These difficulties were attributed to the possibly erroneous position- 
ing of Point 1, the possible misidentification of Point 3, and the uncertainties 
associated with the coordinates of all points that were fixed and employed as 
control for the model. Elevation differences were particularly noted to be 
a potential source of error since the elevation differences among all points 
were relatively small and generally within the predicted error of the lunar 
map (6 meters with 90% probability). A further complicating factor involved 
with the uncertainties in elevation determination and the minimal differences 
was the near coplanearity of the control. As explained by Smith [29] this 
would manifest itself in the triangulation program as an indeterminacy of the 


normal coeffecient matrix. Of possibly worse consequence is Thompson's[31 ] 


18 





expansion of Smith's explanation Which would indicate, if not indeterminacy, 
then an instability of the solution. | 

Triangulations numbered 27 and 29 provided the most consistent 
adjustments that could be extracted from the real data. However, No. 27 
utilized Point 3 as a control point and, as a result, the values involved must 
be suspect. Triangulation No. 29 utilized Points 2, 4, and 5. Since these 
points appeared to be well identified, were in proximity to the camera 
exposure Stations, and had the best known positional Poleatane ip) this trian- 
gulation is accepted as the most valid. Unfortunately, acceptance of trian- 
gulation No. 29 positions the LM at LAT. 3-11-51.4S LONG. 23-23-08.1W. 
This reduces the distance between the LM and Surveyor III to approximately 
163 meters, and redefines the bearing to the LM to about 311° 30’. This 
possible redefinition of the scale and orientation of the system effectively 
distorts the information it produces. Nevertheless, the results do have 
value insofar as the adjustment retains consistency and merely lacks a valid 
scale and orientation. A complete summary of the results and the statistics 
of these adjustments are provided on pages 25 through 36. 

It can be concluded that the possible gross uncertainties of this 
particular set of real data negate any reasonable expectation of significant 
results. However, the feasibility of employing real data with proper control 


seems to be reasonably apparent. 


4.2 Idealized Data_Experiment 

4.2.1 Procedure 

In order to provide a standard by which one might logically anticipate 
the predicted accuracies of a triangulation program utilizing lunar surface 
photography with realistic control, an idealized model was constructed (Fig- 
ure Ill). This model presupposes a reasonable capability of determining the 
relative elevations of Point 3 and the camera exposure stations. with respect to 
Points 1 and 2; and an ability to make an estimate of the xand w rotations of the 


elements of exterior orientation. 
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I, (110, 110, 0) 
II. (116, 110, 0) 





ma (lula eNO} 
Preis, 15256) 
, (118, 132, 2) 
. (102, 160, 8) 
, (160, 300, 12) 


or WwW Fe 


Figure III 
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Additional conditions and parameters are: 


~~ 


1, Points 1 and 2 aligned with local vertical at a fixed distance 


(In this case 6 meters) 


2. Point 3 at a known elevation relative to Points 1 and 2 and 


at a known distance. 


(In this case 5 meters) 


3. A Hasselblad camera (described previously) utilized for at 


least two exposures providing a stereoscopic pair at a 


distance near its 22.4 meter focus. 


The coordinates of Point 1 (assumed to be an LM or other similar 


landing vehicle) are considered fixed and to define the origin of a local cartesian 


coordinate system. Points 1 and 2 define the Z survey axis; and Point 3 is 


then defined to be on a line parallel to the X survey axis. 


The coordinates of 


the control points (1,2, and 3), the photoidentifiable features (4 and 5), and the 


exposure stations become: 


oo fF Cc DO & 


POINT 


(Target) 
(Target) 
(Target) 
(Feature) 
(Feature) 


Exposure Station 1 


Exposure Station 2 


a€ 
113. 00 
113, 00 
118, 00 
102. 00 
160. 00 
110, 00 


116, 00 


» 
132. 00 
132. 00 
132. 00 
160, 00 
300. 00 
110, 00 


110, 00 


And the rotational orientations of the cameras is: 


G 


LOCAL CARTESIAN COORDINATES (meters) 


Z 


. 00 
6, 
oe 
37 

12. 

. 00 


00 
00 
00 
00 


. 00 


CAMERA ORIENTATION ANGLES (DEGREES/RADIANS) 


Ce) 
Exposure Station 1 0.00/0.0000 0.00/0.0000 90.00/1.5708 


© 


Exposure Station 2 0.00/0.0000 0.00/0.0000 90. 00/1.5796 
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The estimated standard errors associated with the various observations are: 


Photo coordinate > Ox = Gy = 0.005 millimeters 
Survey coordinates 
Targets : a, = a, = a, = 0.01 meters 
Heaubes mg yeaeO, — OC, — “Inerers 


Exposure Station : 0, = dg, =°meters 


ao, = 0.01 meters 
oy = 0.08727 radians (5°) 
Oy = Oy = 0.17453 radians (10°) 


4,2.2 Idealized Data Results 


The initial adjustment of the idealized data was a slight modification 
from that which is tabulated. The first triangulation constrained the survey 
coordinates in relation to the relative errors of the lunar map. This was 
assumed to provide standard errors of gd, = Gy, = 3 meters and a,, = 6 
meters. | 

Although the results of the first adjustment produced smaller standard 
errors in the adjusted coordinates of the photoidentifiable features, the realism 
of estimating the X, and Y, of the camera exposure stations on the lunar sur- 
face to that accuracy appeared questionable. On the other hand the estimate of 
elevation differences between the camera stations and Points 1 and 2 toa 
reasonable accuracy seemed practicable. As a result, constraints on X, and 
Y, were removed and that on Z, was strengthened. These results were predict- 


ably good and are provided on pages 37 through 42. 


In an effort to produce results that might be more indicatory of 
those that could be achieved in actual lunar surface operations, the coordinates 
of the surface features were perturbed within the limits of the map accuracy. 
Expectedly, the results were identical. In a subsequent adjustment the 
constraints on the photo-coordinates were relaxed; that is, the weight on 
photo-coordinates was reduced from 40,000 to 10,000 (,,, = 0.019 millimeters 


vice 0.005 millimeters). This caused a significant deviation of the adjusted 
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coordinates of the lunar features from the known positions. In turn, the 
constraints on the rotations of the exterior orientation elements were relaxed, 
and the camera constant was perturbed by an additive 0.050 millimeters. 


The following table provides these results for comparison. ° 


Condition I: Coordinates of Points 1,2, and 3 constrained to 0.01 meter; 
Zo to 0.01 meter; xto 5°: (and w to 10°: photo-coordinates 
to 5 microns; and f = 60.0 millimeters 
II: All of the above except photo-coordinates constrained to 10 microns 
III: Same as II except constraints on % ©, andw removed. 


IV: Same as III except f perturbed (f = 60.050 millimeters) 


CONDITION | ADJUSTED COORDINATES | FEATURE POINT NUMBER 


ee he i. 0c 





4 

I 102.009 | 159. 991 15995 4 

II 102.010 7 }-159.977 1. 996 4 

Ill 102.010 | 159.977 T2396 4 
IV 102.010 | 160.001 7.996 4 
KNOWN 160. 00 300.00 12.00 oD 
I £597 9904) 2992503 o 

IT 199.907 | 299.514 o 

Il 1597908 | 299, 516 oD 
IV 159.908 | 299. 657 o 





It can be seen that the most significant deviation of the adjusted 
coordinates from the known coordinates of the feature occurs as a result of 
relaxing the constraints on the photo-coordinates. This is not unexpected 
Since there is a large weight change involving the elements which provide 
the basic control for the model. The only other significant deviation is 
noted when the focal length of the camera is perturbed and this is apparently 
confined to the y survey coordinate which coincides with the rotated camera 


Z axis. 
2o 





Although the scope of this investigation inhibits specific predictions 
of accuracy, it appears that with proper control on Points 1,2, and 3 and the 
Z, of the exposure stations a calibrated camera is capable of producing : 
positional accuracies of lunar features to several tenths of meters at distances 
of approximately three-hundred meters from the control. The limited number 


of points negates any empiric estimate concerning the relationship between 


positional error and distance from the established control. 
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JOB NUMBER 0 
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RESULTS 
SURVEY COORDINATES 
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5. CONCLUSIONS 

It is apparent from the results of the real data experiment that, a 
general, the potential to improve selenodetic control by the use of lunar 
surface photography exists re a significant degree. Algeanen the specific 
results are considered inconclusive because of the lack of any dependable, 
local control, the experiment has emphasized some of the difficulties 
associated with surface data. Of particular note, is the instability of the 
solution due to the relative coplanearity of the control utilized. This is 
a realistic problem when one considers that the APOLLO landing sites to 
date have been selected in the mare areas where relatively level lunar 
terrain has been a criteria. It is anticipated that this criteria will 
continue to be considered, but perhaps, to a lesser extent as the expe- 
rience in lunar landings is increased. This does, nevertheless, stress 
the need for good vertical control, strongly constrained, to minimize this 
instability. Additionally, the solution has manifested a certain sensitivity 
to the rotations of the camera's elements of exterior orientation. This 
was particularly evident when all elements of the exposure station were 
constrained and Points 2, 4, and 5 exercised total control of the model. 
The resulting adjusted coordinates were realistic only for those points 
and stations within approximately fifty meters of the control points. The 
exposure station for photo AS12-48-7090, the most distant of the exposure 
‘stations, was almost two hundred meters from its estimated position with 
more than twenty times the estimated x rotation. Point 3 could not even 
be plotted on the chart. 

On the other hand, the idealized data and that with perturbations provides 
some indication of the kind of accuracy that may be achieved by a reason- 
able effort to establish a local network to control the adjustment of more 
distant features. Further, one may reasonably imply that an additional 
input of data from overhead photography (properly scaled, if a camera 


lens of different focal length is employed) would provide a material 
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improvement to this adjustment. Yet, no specific predictions can be offered 
because of the paucity of points and the lack of suitable overhead photo- 
graphy and information regarding the lunar conditions (such as surface 
refraction, etc.). However, it is justifiable to assume that the photo- 
grammetric errors associated with the adjusted local eoemates of lunar 
features from surface and overhead photography would not contribute 
materially to the total error substantially attached to any astronomic 


observations. 
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6. RECOMMENDATIONS 


6.1 Concept 


The procedure to be described is a direct application of fundamental 
geodetic and photogrammetric techniques as described in most textbooks 
on the subjects. The unique aspect is that these techniques are applied to 
lunar surface photography supplemented by orbital photography. The 
basic advantage of this proposal is to establish control where selenodetic 
control ought to be established...on the lunar surface. 

Although this control will be limited in coverage, each subsequent 
landing will provide a further expansion of the control net with an ever | 
increasing number of points which can be identified from orbit and to 
which a set of astronomic coordinates originating with the LM may be 
associated. It is theorized that eventually a net of sufficient extent would 
be available to effectively control unmanned, orbital photo missions. The 


following procedure is offered to that end. 


6.2 Presupposition 


The current lunar landing vehicles are capable of obtaining the 
astronomic position of the landing site from stellar observations. It is 
presupposed that this capability will continue and perhaps improve in the 
accuracy of determination as the APOLLO series progresses. It is further 
assumed that an azimuth can be determined to relate any local coordinate 
system to the selenographic system. One method that suggests itself is to 
image a Stellar field on the lunar surface photography related to Universal 
Time through spacecraft time. This might be accomplished by the use of 
a half-circular, neutral density filter for the Hasselblad camera. The top, 
or clear half, would permit sufficient exposure to image the star field while 
the bottom, or tinted half, would inhibit overexposure of the lunar surface. 


Time of exposure could be recorded on a magnetic taped voice circuit. 
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6.3 Equipment 


The following equipment is additional to what is carried on the lunar 


module, and serves only as an example to accomplish the desired 


procedure. 


i. 


A calibrated tape of approximately 6 meters that can be 
hung without interference from an available or added 
projection on the LM. This tape would be targeted at 
each end with an additional target whose position can be 
varied and its reading noted. A second, similar target. 
for exposure station reference is optional. It is visualized 
that they would slide on the tape with friction clamps to 
maintain their position once established. The lower end 
should be weighted and might have some dampening device 


to reduce oscillations. 


Two lightweight, variable height, telescoping tripods. 

a. One targeted tripod would be equipped with a 
small leveling telescope and two calibrated, 
horizontal spirit levels. One glass parallel to 
the telescope optical axis, the other normal to it. 
A plumb bob or optical plumb is necessary. 

b. The second tripod would provide an attachment 
for the Hasselblad camera with similarly 


oriented spirit levels. 


A calibrated tape of convenient length (perhaps up to 20 
meters) with staking rings at each end and a tension 


spring with scale at one end. 


- 6.4 Procedure 


At any specified time during lunar surface excursions, the astronauts 


would carry out the following procedure: 
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I The vertical tape would be hung from the LM. 


Dr Within 20 meters of the LM on reasonably level terrain 
that would include a background with a maximum number 
of discrete features, set up the level telescope in such a 
manner that it is level and its field encompasses a portion 
of the vertical tape. Position one of the adjustable targets 


so that it is centered on the telescope cross hairs. 


36 Lay out the 20 meter tape from the base of the vertical 
tape to a position below the plumb of the level telescope. 
The tape may be staked in position with a predetermined 


amount of tension indicated. 


4. Position the camera at its first station such that its optical 
axis iS perpendicular to the vertical tape, though not 
necessarily in the same plane. Include in its field of view, 
the vertical tape, the targeted level telescope and the desired, 


discrete lunar features. 


oe When the oscillations of the vertical tape are minimal, 

expose the plate and record: 

a. The reading on the 20 meter tape below the 
vertical tape. 

bi All spirit level bubble positions. 

C. The reading on the 20 meter tape below the 
level telescope plumb. 

d. The readings of the variable target(s) on the 
vertical tape. (All readings could be voice 
recorded on tape.) 


For subsequent exposures, it would only be necessary to reposition 
the camera to obtain a stereoscopic pair, possibly readjust the optional 
variable target (if used) and to record the readings already mentioned, 


(See igure ITI). A” 
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The resulting models would be similar to the idealized model described. 
The vertical tape would provide the scale to an estimated accuracy of a 
few millimeters and define the local vertical as the local Z survey axis; 
its lower target would establish the origin of the local coordinate system. 
The targeted level telescope, the position of the variable target on the 
vertical tape, and the measured distance to an estimated accuracy of .01 
meter, could be computationally corrected to define a line parallel to the 
local X survey axis. The local Y survey axis would then be defined. 
Reduction of the recorded readings would give the survey coordinates of 
the level telescope target and the Z survey coordinate of the camera 
stations to an estimated accuracy of 0.01 meters. The X and W 
rotations would be estimated to be near zero based on the camera 
leveling results. The © rotation would be estimated by its relation to 
the defined YZ survey plane. Approximate positions of the lunar surface 
features can be scaled from a convenient lunar map. 

After preprocessing the necessary information and providing the photo 
coordinates from COMCORDON to the BLOCK TRIANGULATION PROGRAM, 
the resulting adjustment would photogrammetrically relate all discretely 
imaged lunar surface features to the position of the LM in a local cartesian 
coordinates system. 

Extending this with an azimuth and the astronomic position of the LM, 
this adjustment, with a simple coordinate system transformation program, 
would provide selenographic coordinates and relative elevations of lunar 
features that could be related to current and future orbital photography. 

It is acknowledged that the foregoing method is neither the most 
simplified nor the most sophisticated that could be employed. However, 
it does serve to emphasize the fundamental requirements of the system; 
that is, the establishment of adequate scale and orientation and the 
application of sufficient constraints to obviate coplanearity of the model. 

It is, therefore, suggested that the previous procedure, or one 
fulfilling the same basic criteria, be considered for adoption. It is 
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firmly believed that its implementation would be the beginning of an 


improved selenodetic control network. 
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APPENDIX I 


COMCORDCON 


The purpose of COMCORDCON, an acronym for Computor Coor- 
dinate Conversion Program, is, as its name implies, to convert stereoscopic 
comparator observations of object space points on a photographic plate to 
photographic coordinates in a photo-coordinate system. The resulting measure- 
ments, sequentially corrected for film / emulsion shrinkage and radial lens 
distortion, are then in the proper preprocessed form as partial input to a 
simultaneous block triangulation program to be described. 

To achieve this objective COMCORDCON utilizes in this case negatives 
or diapositive plates from a calibrated HASSELBLAD camera equipped with a 
focal plane reseau grid and nominal 60.0 millimeter focal length lens; the 
Zeiss, PSK, Precision Stereocomparator augmented by an ancillary IBM 
026 card punch; and the OSU, IBM 360/75 computer using a FORTRAN 
IV G, Level 18 Compiler. Comparable systems may be substituted with 
necessary modification of data format to provide compatibility. 

Input consists of the number and interval spacing of the reseau marks 
on the focal plane grid; the type and identification of the plates; lens distor- 
tion values; the coordinates of the center crosses, the object space points, 
and four neighboring reseau marks (preferably surrounding the object space 
point); and other job keeping information. This information is then provided 
to the SUBROUTINES TRANCP and RADIS 1. 

TRANCP performs an affine transformation by solving the trans- 
formation parameters for each set of four observed reseau marks associated 
with an image point by a least squares adjustment. These transformation 
parameters are then employed to obtain the photo coordinates of the image 
point, corrected for film shrinkage, in the reseau coordinate system of the 


photograph. The mathematical formulation for this transformation is, for 


D9 





example: 
Assumptions 1. Reseau interval is exactly 10 mm. 
2. No correlation between observations of unique points and no 
correlation between measurements of x and y coordinates 


in an image point observation. 


y! 






RESEAU 
_-® 
. 


(FOR NEGATIVES) 


COMPARATOR 





=A 0h + AX ot A oey’ 
Va Ose Bt xX’ 4.82 yy" 
F, = X-A0-Al1l X‘'-A2Y‘'=0 
Fy, = Y-BO-B1X‘’-B2yY‘'=0 
ae ie Ee 
~] -x} -y! F 
px- |2*| - ei EX=|F =; * 
0 A ~l -X, -Y3 Le 
a Se menyy OBS F 
4 x4 
L. 
OF 
py. [22 E Y = F| 
08s 
vo 1c — bb: ix = BEX; by = BEY 
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PARX = -C'UX 
PARY = -C' UY 
A= A+ PARX s 
B = B+ PARY 


The image photo coordinates corrected for film shrinkage in the reseau 


coordinate system are then: 


x’ = Pom ens Omi en( Ge ) 


A2Bl - Al B2 [Where X and Y are 
the average of f 
,_ BI(X-A0) - Al(Y¥ -Bo) ee ames 
Vee observations. ] 
A2Bl - Al B2 


Assumption No. 1 is not necessary since the reseau marks can be measured 
individually and their true coordinates read into the program [17]. The 
standard error of their position can be determined and propagated through the 


complete system if desired. 


Assumption No. 2, though not entirely true, can be minimized by careful 


attention to proper observation techniques. 
Both assumptions are made for the sake of simplicity. 


SUBROUTINE RADIS 1 employs the current image coordinates and applies a 
suitable correction based on camera calibration data for radial distortion as a 
function of the radial distance from the center cross. It assumes a symmetric 
radial distortion and that the center cross is coincident with the principle 
point on the plate [17]. Again those assumptions are made for the sake of sim- 
plicity. It would be entirely reasonable to integrate into RADIS 1 a correction 
for the deviation of the center cross from the principle point, a model for 
non-symmetric radial distortion or a thin prism model for correction of 


decentering and radial distortion [11] [10]. 
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COMCORDCOM 


LIST OF PERTINENT VARIABLES AND INPUT DATA 


Card No. 


Variable 


CR 
A,B 
INFO (1) 
INFO (2) 
Cont) 


CC (3) 


Ce (>) 


CC (4) 


COND 


TY PE 


INFO (4) 
INFO (5) 
CCO (1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(12) 


Purpose 


True Reseau Coordinates 
Transformation Parameters 
J ob IdentificationNumber 
Reseau Interval (in microns) 


X coordinate of left plate 
center cross (mm) 


Y coordinate of left plate 
center cross (mm) 


X coordinate of right plate 
center cross (mm) 


Y coordinate of right plate 
center cross (mm) 


Defines the plate (s) being 
used, (RIGHT, LEFT, or BOTH) 


Defines the type of plate, 
(POSITIVE OR NEGATIVE) 


Left plate photo-identification number 

Right plate photo-identification number 

Point Number 

Point X Coordinate, Left plate 

Point Y Coordinate, Left plate 

Point X Coordinate, Right plate 

Point Y Coordinate, Right plate 

Blank 

Point Number 

X Coordinate, Upper Left Reseau(#1) Left plate 
Y Coordinate, Upper Left Reseau(#1) Left plate 
X Coordinate, Upper Left Reseau(#1) Right plate 
Y Coordinate, Upper Left Reseau(#1) Right plate 


Blank 
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Card No. 
4 Repeat Observations of Card #3 but Observing Upper Right Reseau (#2) 
5 Repeat Observations of Card #3 but Observing Lower Right Reseau (#3) 
6 Repeat Observations of Card #3 but Observing Lower Left Reseau (#4) - 


The set of cards (3,4,5, 6) is repetitive by N number points. 
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COMCORDCON 





GOut , IMPLICIT weayes lA-H,(=7) ; 
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_—- aE se ee 
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221) #7 INtGOS)=NPCINIA] 7 ar eee eh O09 ho. pees 
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APPENDIX I 


BLOCK TRIANGULATION PROGRAM 


The Block Triangulation Program performs a spatial resection and 
orientation by a simultaneous least squares adjustment treating photo and 
survey coordinates and the elements of exterior orientation as observed 
quantities. As significant input it utilizes values of the camera constant 
(focal length); corrected photo coordinates output from COMORDCON; 
observed values of exterior orientation equated to the initial approximations; 
and various job keeping information such as number of photographs, number 
of unique points, photograph number, etc. Additionally, variance-covariance 
matrices are estimated and input to the program for the photo coordinates, 
survey coordinates, and the elements of exterior orientation. 

The simplifying assumption that there is no correlation between 
photo coordinates, survey coordinates, and the elements of exterior orienta- 
tion is employed. It is further assumed that no significant correlation exists 
between the individual survey coordinates, the individual elements of exterior 
orientation and individual photographs. 

The output of this program provides the minimum variance solution 
for the adjusted values of the elements of exterior orientation and the survey 
coordinates, 

This solution is effected through the following mathematical form- 
ulation: 


The observation equations are: 


6 8 8 8 
V+BA+BATtT€E = 0 
6 6 6 

2 ee LS +€ = Q 
8 8 8 
V-A +€ = 0 
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In matrix form 


WW oe 
KW 
CS o 
mn 


V 

@ 

Vi + J-I OF} AY + 
gs 

V 0-I 


~ a™® @ 


+ €= 0 


eel 
D1 


vee 
Then by exercising the method of LAGRANGE for realizing 


the minimum variance solution we have: 


thus 
(B'WB) A + (BWE) = 0 EQUATION A 
where 
e 6 0 Xy ae 
B — B 8 . B= B B, 0 2B = 0 (Xo, Yos Zo, Hos Poy Wo 
B ; cee ats aY 
Pen O( Xo 9 Yos Zo» Hos Pos Ww o) 
i B sae . = Vaz 
B = B 8 > B, a : 
Oo 'B oy 
ae a(x, Y,Z 
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W,, 0 


: w °0, WwW, 0 
ae ° \ = ° _ 
W W Me We ; Ws 0 Wyy 


W 0 


0 Wuowo 


W xx 0 
W 
0 es 


ar . |e aaa a 
A = A ; A = De ; 5 A = . 5 6, = bY; 


A € 
5s 53 
b e ° 643 
65 3| bon 
& 
€ ese 5 Ky = Ry 
- e 6 = 3 
te =; y\5 ett; 
i ° = G = Wy 
€ Yn — Ya ~ 


O 
8 
I 
F 


oO 


i 
J 
8 hg 


oO 
Y;- Y 
Zin 


8 


WY 


© 
Zan- Zan 


—_— = a Ul ee 6 64 
BWe = U = . Ui = Blwe -We 
U2 8 8 8 

| U2 = B'We -We 


with m = the number of photos and n = the number of unique points, 
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Expanding EQUATION A provides; 


T oli/B Bille BT-I 0 0 . 

. We, So weal tcl . 

; eealiale 0 : ; Wiss € =) 

B’ o-I{] 0 0 =] |i B' 0 -I 0 : 

and after multiplying out, 

8 8 8 8 s 8 8 8 

Bev B+ Wi BWR i ‘i B’WeE-WeE 

as ”6hCCOC~<“ Ct aves = 0 EQUATION B 
B'WB iBT WBtwW iiiA B'Wwe-we 


By inspection one can ascertain that Bt W B+W is a full six by six matrix with 
one photo, but for m photos, a block diagonal matrix of 6 m by 6 m consisting of 
six by six submatrices. Analogously, BT W B and BT W B will be 6 by 3 n and 

4 n by 6 posecelyely for one photo but 6 m by 3 n and 3 n by 6 m for m photos. 
BTW B + W will naturally be 3 n by 3 n for one or m photos. However, 

because of the assumption that there was no correlation in W and W, BT W B te W 
is a block diagonal matrix of three by three submatrices. These facts concern- 


ing this system of normal equations lend EQUATION B to a simplified method 


of inversion by partitioning [13]. 


8 8 8 8 8 

BIWB+ W | B'WB A] B 
i | eo abe (= 

B'WB | BTWB+W Culucl) 


where N is of the rth order and A and D are of order 6 m and 3 n respectively; 


r=6m + 8n 
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therefore > oie BM 


It can then be verified that; 


Ko (A= B De): 
he te Bebe 
ies Derk 
Ne= peepee i. 
. = ! ; 
Coa (A -BD“‘c)*;} KBD? 
PN aa a (pee tle 
-_D'C K : D'-D°C L 
6 6 6 8 8 8 Cc 8 
and where N = B!' WB; N = B' WB; N= B Wes 
6 8 ¢ 8 8 c ma Cc 8 8 
(N + W)- N(N+W)? NT | -K N(N+ W)? 
N? ae ae /_ sO + en eee 
8 8 ; Cc | 8 8 r 8 8 : oe 
-(N+ WwW)? NK ) (N+ W)° -(N+W)* NUL 


As a result only one 6 m by 6 m matrix, andn three by three matrices must 
be inverted in the partitioned matrix to effect the inversion of N, an rth order 
square matrix (r = 6m + 3n), for m photos and n unique points. In solving 
the normal equations the solution then to the alterations to current approxi- 
mations is: A = -N7*U and the adjusted solution vector is computed as; 
neel= x + A 
with subsequent iteration converging to a suitable solution. 
Kage = Aa, + A, 
The BLOCK TRIANGULATION PROGRAM is capable of processing up to 
twelve photographs, 36 unique points, and a selected number of iterations 


with but minor modifications. Additionally, it provides the standard error of 
(i 





unit weight, and the residuals, standard error, and the variance-covariance 
matrices of adjusted parameters. 

Other than standard library subroutines this program mnie 
SUBROUTINES MATINV and COFEI. MATINV is merely a double precision; 
matrix inversion routine employing a bordering technique; COFEI computes 
the partial derivatives of the projective equations with respect to the elements 
of exterior orientation. The main program provides for updating photograph 


and point information for use on successive passes of the iterative process. 


The terms in the preceding discussion are defined as: 
V = residuals on observations 
B = partial derivatives of the projective equations 
4,5 = alterations to current estimate of the unknowns 
¢ = functional relationship between the observations and 
unknowns evaluated at the current estimates 
X= Lagrange multiplier 
W = weight matrix 
x, y= photo coordinates 


X, Y, Z= survey coordinates 


Xo, Y5; Zo > 
x ,P,W 


= elements of exterior orientation 


Superscripts 

©. refers to elements of exterior orientation 
; refers to survey system 

none; refers to photo system 

; refers to observed values 


©; refers to computed values using approximations to unknowns 


- 
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BLOCK TRIANGULATION 
LIST OF PERTINENT VARIABLES AND INPUT DATA 


CARD NO. VARIABLE PURPOSE 
A Elements of N and N* 
U1, U2 Constant vector 
DE Alterations to elements of exterior orientation 
DS Alterations to survey coordinates 
E Residuals on x, y photo coordinates 


sin(xX, O,W) COS( xX, Mw) 


1 CC Camera constant (f in mm) 
N Number of photographs 
2 WM A scalar multiplier of a 2 by 2 identity matrix (I) 


to provide an estimated variance-covariance matrix 
for observations on photo-point coordinates 


IP Total number of unique points 

3 Photo(1, 1) Number of points in photograph 

4,5,6,7,8,9 Photo (1, 2-7) Current approximation and observed values 
(X0, EE] (equated) of the elements of exterior orienta- 


tion (m, and radians) 


10,11, 12,13,14,15 Photo(114-19) Estimated variance-covariance matrix for 
exterior orientation 


16 Photo(1, 1,1-2) x,y photo coordinates (COMCORDCON) 

7 Photo(1, 1, 3) Point number 

18,19,20 Point(1,1,17-22) Current approximation and observed values 
(X, ES) (equated) of survey coordinates (m) 

fl 22,23 Point(1, 1 8-10) Estimated variance-covariance matrix for 


survey coordinates 


CARDS 3 through 15 repeat by m number of photographs and 16 through 23 


repeat by n number of points in each photograph. 
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MAI 


BLUCK IRIANGULAFION 








Cc LECK BLELBEEULLEBbIGE oes 
- Se ia en as ee ee Ee ae st ee 
ooul IMPLICIT REAL ENC A-H,O-2) 
meee —leti(i‘“‘;‘;« ENS ALL 018 1H) AL 20619,125) -A221156135) ' 
YOU 3 UVITEGSLOW BEATCL belo) AZZELTS TU OECLY) USIL5) | 
~~~ QO0U4 a CLACNSi ON ULL1G) ,U2E19),AnO 18) PHOTO 3,49),P0INT(3,5,22) ' 
C005 UdhewSiON XO(6),-CF(0).nWE( 6,6) 
“Q0U0u “GURUS TOs DATAC2 Dewl et e2 eXC Fd ecSl 3) owSl 3, 3) ° ry 
vous DEMLNSCUN LEAR AATLMMN (Ge 2) U2 a9 RIGO) SMONIH(L2),1TITLELZ) 
~ 9008 DIMENSION KLELOO) sKLL (100) 
0009 INICOLR PP,P,1T,KEY(35)/35*07 
—oo10 ~~ "~~ DATA BMUNDHELIS 9 SA. Os eMOUTH( 2) "FEB. 9 / MONTHI 3) /*MARS S/S MONTHI(G)D/? 
PAPR ES SMUNTHCDISEMAY "/,MUNTHIOID/S* SUNE*S PMONTHITIS*SULY'/eMONTHIS) 
——  ltlttttt—S S/S AUG A eHONTHLIDIS*SOPT YS MONTH LOIS ECCT. YS pMONTHILIDS *NGVe 8S ,MONT 20” 
YHE LZ TCLC. M7 
— GOL 7 DALLA RM ARNG KIGSYEMY, MND, OSES 
0012 QOQATA TELK/° */ 
0013 oo TOUT=0 
0014 UF=C.0 
im OO15 Vn en re eh Sa ht i ks oo Cp 
OULG ICYCLE=0 
our €6=6—0Shftéi‘(i‘éié™O...  ©6©63. US : os . 
OOI8 REAL( 52919) (TLILE, JOBNUM 
OULD = 7 ReAddlo,2s9) CCov, 1? 
cCoz0 REAU(5»2799) wm ile ee 
—0C7, = =———COM I FGA AT CF LO53,275)—O an’ 
OU22 DO 2CU 1=1,N 
~ QOO23 —™  " READL5,603) PHOTOII, 1) ~ 
OUZS P=PHOTUC(,1) 
Go7> °° °° ©43xG 201 Jsa2,7 " : 
0026 KLEAU(9~5900) PHOTOLIL,J3),PHOIO(1,3+6) 
Seow ss 8=©=©=©—CU OT PHOTOCI, J+6)2PHGTOLL.J) COCO ~ 
0028 L2=13 i ; 
“cuz?7 - °° °° «4400 202 J=1,6 
0030 “ tleL2el _ . . 
UO3slL 7 ORD EL 246 . ae ~ 
00 42 202 REAL(5,501) (PHUTUCI,K),K=L1,L2) 
~ GO33 DU 200g athe eee 
0034 READ(5.,9G62) POCITIL,S»,LI,POINT(3,5,2) _ 
' 0035 ~~" ~~ "REAGI5,663) POINT(I,J,3) as _ or 7 - 
00 36 PUINT CL etde4) en 
Seocs,  4£4POCINICT,J,5)=6.0 
0638 POINTII,J,6)=0.0 a 
soe  ..  4.£POLNT(L,J,/)24¥ #  &@&@&+«”}”=—"” +  @©~©6C.WhCODrD)KmKmW©3©€~ : 
0040 GO 204 K=IT7,ly 
0041 ~~ © *  REAG( 9,969) PUINTOCI,JeK),POINT(I,JeKe3) 
0042 204 PUINT( 1, JeK*+3)=S=PUINT(I, 5K) 
RT re 
0044 . OU 265 K=1,3 
~~ 0u45 ileal : 
0046 L2=L2+3 
—— 0047 205 REACI5, 504) (POINTOITZ JIL Lely, ta) Sia eracy evap ee 
0048 200 CONTINUE 
C LLEAK NURHALS 
~ O04) KT=3* TP 
0050 J1=608N ale 
mous, 6 §EYCLESICYCLE¢] ~~ ~~” ‘ 
0052 DI) 5 T=1,1P rs 
f O53 ~ ~~ 5 "KA EVEL) =0 OR ene Se pe ee oe a eT Se oo 
0U54 LUC r2 1 =1s 31 
SoCs £  UOt0{y=O.0 
0056 VU ¢ J=l,J1 
Se GOSt) °° 4g ALLC1,J3)=0.0- 
0958 OQ 2 J=1,K1 an 
ee a eee le) 0.050 = ne eee 
0060 VU 3 1=1,K1 
S O06) U2) =O. 
0062 0G 3 J=1,K1 ‘ 


Meo los6 oe oes A22tl,J)20.0° °° ~—sCCOC—C—C—COC 
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wwe oe ee 





Cease HoTGcOALTA CCC 
G 
meeooi.  . . . £4§'Pcsl ae 
C0605 LO 100 J=lotd 2 eee 
“06 °  &°3tCC=)3 ie aes ae 
00u7 Ke=PHOTCT Jel) ee ee 
woo: uv 210 [L=2,7 —a lll ll 
GVO9 AGU LTI-LI=PHNTEII LD) 
~ 0070 EEC T(-lbePndf00CI,TT+6) 
0071 VO ZIG JJ21e6 ig a Cee ee oe 
= soo72 = ~~ LOCFLOC +) ic °c pares 
0073 210 wEELI- yds) =rtCTUI LOC) 7 ke 
70074 °" ~ WRITE Lo.GO1) XU aera ; i : 
cols walleloruIt) KES 
eevee.) OOL FOAMAT(6E15.5,/(6E15.5}) °° & &«@ & &  . @©@]|+6»9"D’”OC 
00/7 = OS liv eKl 4) ) ae a ; 
Severs) j. £REZVENGINIXOIS)) ~~ ~ i 
0079 RO3PVZUSTNEXULG)) 
“co3.0 @#@#&«;* # # ;2§ K(4)2DC0S(40(4)) - 5 
0081 Ri>o)euCUSI%uUI95)) ae ae. 
“—egoue  . .. £Kt6ENVCOSCXOTG)) ie a 
C 
—< Cc — ADDING TERMS OLE YO CONSTRAINTS ON EXTERIOR GRIENTATION ELEMENTS 
C e 
0UusS3 ~ po eel is 1) a ~ " ‘ 
0054 vO 10 [1=1,6 


~~" C085 TF (aFUCL LEY NEL OO) 


a S 


OF2vUF+1.0~ 











© om ee em eee = 


—~ ema o =-—mecce -—— eae - 


ee ee 


—— eee ee 


2 6 eee 


0086 VUFFVDF*enEC( List Lye(CElLlLiy—XOCT1)) 442 
—veat °° 4 NG 1) JJ2l eo 
0064 VUIILe LL FULILe DD -wWECTL JD FeLEEC JI -H XO) ) 
~ OO8D TT PO ALLILe LE Le JJ AALLOLeL I eh ed) eWEC LI, S) 
: Cc 
— +. C CALL POINT DATA 
C 
- 0030 ~~" ~~ vO (GO L=t,K 
0021 DATALLI=PULULI Ie fel) 
oo32 =" ~~ . DATAI2) =PUINE (Se4e2) - 
0033 P=PCINT(J,1, 3) 
~~ 0094- Loc=s3 = el 
COX UG 215 I[k=1,2 
~~ 9096 ~~ VI 215 sJ=l 22 - ; - Fe ia 
Co9T Loc=c_cCe] 
00%. “" 215 wUll,JJ)=POINT(J,1,L0c) . 
C099 UU Z2lo Iz 
~ OLo0 00. 7) Gael) = eee COO 
0101 LOC FLUC #1 
Gee” —6060—CCt”t~—~—“—;~CC*idCK OW SCT LAI EPOINT (Je 1,t0C) 
C103 DO 217 [1217.19 
0104 X(EL-LoP=POINT (SeT, ST) 
0105 217 ESTII-16) sPULNT (dele 1343) 
~ O1L00 wWRITELosoO0l) KX 
O1L07 wRITELGsu01) ES 
- O63 ~~ —" DF=DF 42.0 
0109 LL=3e[r=-1) 
o : -¢-- 
G CALL PARTIALS 
ee Cnn ee 
0110 CALL CCFELIXR,B8sXO,CC,XC, YC) 
z —— : 
C COMPUITF NORMALS 
Ct °° °° °° °° (Tesi=l ae ee ase 
U1l12 [Fixcy(P).GT. 0) GO TO 12 
~ 0113 KrYy (2)? 
C114 LTESTtT=0 
Gui 42 00 > [l=1s6 I ee ena 
0116 09 19 JJ=1,2 
"O17 & 3&7 —COOtTOCO:COP mB EMP LLL 2JI) 20.0 - a aie coahs a7. ve 
0116 VO 15 KKz1¢2” 
7 O11l9°~——— TSS ATERPEIS Le IJ) = USS) NU Ue 
0120 CIEL) =UATALI1L)—xXC 
~ O1L200 £02 2D FOATAL2I HVC ce 
0122 _VOFSVUOFeW( Ls LIFECLISFZ Wl 2, 2)9E(2) 482 


_——— +. = 
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Vids 
0124 
O125 
O26 
Gre? . 


16 


0128 ele 


cme se eg 


UlL29 
0 °° ~~ 
_ O13i 


we ee ee ne oe ew es eee 


MAA 


“SMe LS - bel 


~ ENE = 7 49 


™ "8" CumPuft aAa2l 


COMPUTE ALL wt et aa a 
VU lu KK21,2 

Vibe PL sUblbelL LL) *¢ATE APC LI, KA) CEL KK) 

vl) Ef Jbdeb et 

LO Lf AK=1,2 


ALLCLOL Es Ct SS) EALLOLE CT U4 Js) tATEMP CLL KAKI OB (KKe SJ) : 


COMPUTE 412 
UU 19 Ll=le6 


J0- 19 ae 1,2 





ee 





So oe ee TT RE 








CUMPUTE A22 


~ 1s 3 DG=23 1f37,9 


0134 
>. 
0136 _ eas 


GLa eS 
C138 


~~Q1359 Py a 


0140 
~~ QO141— 
2 


0143 26 
0144 


—-~—“C~SrsSs:C‘<CS*~*<;273 3 3S oe + eee ° oe 
e 


Cc 


Se Ce 


ce ee Se ee 


~" y2th)su2l4) +AVEMPELL, KK) EC KK) 


ro 


vi 23 SJ=1e2 
ATevr(ll-6,5J3)=0.0 
UU 23 KK=1.? 
ATEMPCL(-GeJSJIDFATOCMP(L 1-6, 35) ¢Of KK LL IFW (KK IS) | 
DU 2¢ (t=1, 3 


- a ee ee eee a ee 








UD 26 KK21 2,2 


VU 20 JJ=749 
A22(hyLLeSI—G) =A22Z(N LL +JSI—~“6) CATEMP( LL KK) *BIKK, one 
DFCLTtST. EO.L) SE TU 27 


me meee a 





ee =e 


ADDING TEAMS DUE TO CONSTKALNT ON SURVEY CGOROINATES _ 


VDF SVDF eaS CLL stl) eCESELID<XELL) )ee2 - 

LFQnSUUL SEL) .NELO.O) ODFSEDF4l.O0 ~ i, 

OU 28 JJ=1,3 

Rpt TIT=AC2(hsLL eS) ewSC1ll,JJ) °° °&#°«#»°©€©—C 

UZIN) FU2ZUM—wSCLIL JI) FC ESI III HX (ID) 
CONTINUES - Oe ia 
CUNTINUE 
LF(1GUf.cOQ.1) GO IG 300 


ee a ee te 8 re ee ee ee 


— ee ee ee ee 


a.) 0CCt~:*C«SGRPEWNG LIT 


& 
C 


—8 eee eee ee ema se eg ee 


0154 35 
°F °@&°&#7;7&7&7}”}8§§€—CC 


nr ene me 


~~ CACC MATINVOAZ2, KL) 


AL2SUT pe SIFALZC Ts JIMALIUL KI BL 20 (Kd) 


UO 35 121,Kl patie ee a 
vd $% ; dete l 

AzZ2tt{,J)-0.0 

HW Wy REE eKI . 
A2Z2TUL,JI=AA22T (1, JIFAZ2S]UTKIFALZAC Sy K) i . 


OU 30 l=leJlL ia ee 


DO 36 J=leJl 

OG 36 aal,£4l 

ALL (Le JD=ALLUL  JI-ALZU 1, KI #A22T (Ke J) 
CALL MATINVOUATT, J} 


-_— -_ —~ = —« ee -e - . . - 





FORMING ol2 
VU 37 [=1, 51 - 
DO 37 J=leKl 

vl2ft(t,J)=0.6 
VID 3¢f K=1.Kl1 ; 

itt eheel (Le JItAL2(1.KICAZZ(K« 3} ° "= = &~77tCSsi—iai ST /'' dee: 
DU 38 i=leJl 

DO Ja J=l,k} 

Al2(l,59) 70.0 

Ub ts K=is J 1. 


> Saw ee ee a im ee cme ee ee eee oe = 


e 
mere «- - -_-= ——e - wee @+ © meee oe _— es - - - ——_e.: ~- 


Se ae ee ce oe + - tee tere = —e se we e- eso -= 
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C ~ FORMING 822 ° 











O174 OL 39 1=1,k1 
~O17> ~ YU 39 Js1,Kl ae a 
0176 DO 39 k=1,J1 : ; 
of : 39 ASSL SV=AZ20LJI-AZ2T OL, K) FALZ(K J) . - - 
~ "GO" SULVING WORMALS 7 . eee a - 
G ‘ 
Boris ~ remem GM Tea ei ne ee a 
0179 LVE(1)=0.0 
0180 Di) 4] J=1,J1 
O181 ° 4L DELL =UECLI-ALLOL,g) eUl gy) 
ace  — VO 4G J=loAK} 
0183 40 VEC IISUVECILI-AL2(¢1,5) 0209) 
~ 6164 es ens ee aa T= 1s a ae 
0185 OS(1)=20.0u 
Olo&G = = — > 66763 S=14K1 eae nae a > 
G187 G3 USC1V20801)-A2201.5) 80209) 
6168 = TG tee he Vel as ll 
eres | gees SINUS 0a SS a 
a _. © APPLY ALTERATICNS 
C1490 i vu bbs L=L,N Ee Ta a ? 
Olvt RePHIVGTL, 1) a. Fe ae ee 
“Ghue ie a OOOO 
Oly3 PHETOL Ls fl. =P hOTOCL, TZ) eve (6st l=-L)elZ—1) mea... 2 
0194 LOl CUATIALE 
OL1YS VU LLS J21,K a _ 
ClY0 s P=PUL abt leJs3) 7 
OLY7 DO bts 12217,19 — 
“O19 PULUTC Ls Se EL ZIY=PUINT OTe J, IZ) 40S( 390 P=1) 12-16) 
0199 Ll> CUNTINUE SN ais eet. 2: 
a) lt ”:~—<“i;t*sSCSS SM OV CUELUIWG) GO 1G 4  ~ TO CC 
6201 CF=0.0 Pare Sees... 
0202 ~~~" * vut=C.0 rere 2 06 6 yas 
0203 fubl=. 
0204 SiGe ea oa 
0205 3CO CALL MATINV(ALL,J1) ; secon as at = 
0206 — — SRE RViNe 2s ml llC!™C~CS:*C*~<“<‘<i<“i<i‘“a‘ SOS 
020/ DF=LF-Jl-K1 eed - 
0208 TTT oat EV DFR UF as SGM YS "  aviitec 
0209 UNSTUR=OSQORT( VAR) 
“0210 “CALL LUALLMCLYEAR, TMONLH, LOAY,ILTIMED 
G211 Of 346 Lal,J} ee ee 
wiz © OU 30H J=1,J1 . 6 ce =" > ae 
0213 306 ALL(I,JI=ALL(L,S)*VAR le. 2. re oe 
wa.  - * v'} S07 L=1l,K1 a aie ie eee ee 2 er 
U215 Diu 3C7 J=1,K) 
(0216 SUT A221 J V=A2201,5) VAR 
0217 WOF EDF a eet Sl i 
C2138 i UU 301 131.24 ao. | ee 
0219 LFCLIIME.LT. 360000) GU TO 302 ae ee = _. 
6220 — TT) BGt eq hise=36CO00" § ; 
0221 302 LUD 303 J=1,5" 
"0222 ~TFULEIME.LT.6u00) Gu TU 304 
0223 303 ITIMEZISTIKE-0tCO en re ee 
0224 — 304 LIML=FLOATIIII¢ME) 7100.0.  — 
0225 ° t=l-1 | 7 7 me : 
0246 ~~ JeJel , ; 
O22? ARITE(G. 3000) [TITLE » JOBNUM, LOGY» MONTHILMONTH) » LYEARs Lode TIME SNe NO 
na: SF yUNSTOR _ a 5 : 
0228 wRETCO1, 4003) we : ; 
Zee  ° °& & °° °° arnitile, 3005) 
0230 daXilelo,3s001) De tnele ce Ra rn tT te ae eR Fre tn ee eee 
O2at> = = &4+»§ $350. Islan Ce a 
0232 16=o%(l-1) a ee 
02 33 bi 4058121.6°5 + °¢@¢@¢&° °° °&©«& © 
0234 XUCSISUSCRTCALL 1665, 16%5)) ee a ee. 
0235 0 °— SD EEN I SPHUTOC OL eve 7I-PHULUI IL »JSel) 
G236 WRITE(6, 3002) Le(PnCTolled+l) sJeleb)eXO,EEs (PHOTOIL J) Jel ee %9e7) ree aee 


~~ 
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0247 
U2 4 
0239 
0240 
0241 
0242 
0243 
024% 
0245 
0246 
U247 
0248 
D249 
0250 
0251 
0252 


0253 


0254 
0255 
0256 
0257 
0258 
0259 
Q260 
0201 
0262 
0263 
0264 
02605 


D266 
D267 
D208 
0269 
0270 
0271 


0272 
0273 


= 


WRITE(6H,3007) CLALICT6¢K, T64J3), 521,60) -K21,6) 
wWKI bE CO, 5009) 
R(LIZFOSIN(PHUTOCI,5)) 
RIZVFNSINGCPHOTUCE,6)) 
R(JII=OUSEN(PHCTOCTI,7)) 
K(4)=0COS(PHUTOCI,5)) 
RIS)=ODCOS(PHOTO(I.6)) 
R{6)=UCOS(PHUTOC(I,.7)) 
K=PHUTO((,1) 
OO 3725 Jul,K 
P=POINT(I,J3,3) 
KK23¢(P-1) 
OX=PUINICE,J-L7I-PHUIO(L,2) 
DY=PUINE (teste fBI-PHOTOC LT» 5) 
OZ=POINT (e301 9) -PHUTO(1»4) 
XTEVAR(S) FRI GI EDVO(R(GIORILIAR( 3 4R(2)4R14))40Z2¢(R(3) RI 1L)-R16) 48 
(2) *R(4)) 
YT=-DKER(OIFROILICOVE(RIGIFP(S)-“RIZIFRIQIFRILI IV #O2ZE(RI3Z)ER(4)4R(6) 
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